Time evolution of a two-atom dressed entangled state in a cavity 
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We study the time evolution of superposition of product states of two dressed atoms in a spherical 
cavity in the extreme situations of an arbitrarily large cavity (free space) and of a small one. In the 
large-cavity case, the system dissipates, whereas, for the small finite cavity, the system evolves in an 
oscillating way and never completely decays. We also compute the von Neumann entropy for such 
a system, a measurement of the degree of entanglement of the two atoms, as the superposed state 
evolves in time. We find that this entropy does not depend on time, nor on the size of the cavity. 



I. INTRODUCTION 



In the quantum mechanical description of multipartite systems, with Hilbert spaces given by direct products of 
individual part spaces, the superposition principle leads naturally to entangled states which can not be written as 
single products of states of the constituent parts; non interacting subsystems can thus share entangled states that hold 
quantum correlations. Such quantum entanglement carries nonlocal features which can be analyzed by comparison 
with classical correlations p], Q • 

Entanglement is a quantum mechanical resource that plays a crucial role in implementing teleportation of quan- 
tum states and in several applications of quantum computation and quantum information 0, 0, IE]- Quantifying 
entanglement then becomes an important issue which has been addressed in the literature from a variety of view- 
points [E, 0, H, M, EE, EH, EH- For bipartite systems, the measurement of entanglement is well established, the von 
Neumann entropy of the reduced density matrix providing the simplest measure of the degree of entanglement of a 
given state. In this way, maximum entangled states has been constructed for both boson (T3 , Il4j and fermion (l5j 
bipartite systems. 

In a recently reported experiment [la ] , it is proven the existence of deterministic entanglement of separated oscil- 
lators, consisting of the vibrational states of two pairs of atomic ions in different locations. They also demonstrate 
entanglement of the internal states of an atomic ion with a distant mechanical oscillator. The authors claim that 
such experiments may lead to the generation of entangled states of mechanical oscillators in a larger scale, in such a 
way as to provide tests for nonlocality in mesoscopic systems. They also claim that these experiments could be used 
to control quantum information processing based on trapped atomic ions. Previously, an experiment was performed 
in [13], using ultraviolet lasers to entangle two pairs of beryllium ions in an electromagnetic trap. These authors 
also cross-entangled the entangled pairs, that is, entangled each member of the first pair with its correspondent in 
the second pair. Then the first pair of ions was measured, and the results were used as an indication of whether the 
unmeasured second pair was entangled. 

In the present paper we study the time evolution of an entangled two-atom state, in the presence of a force field. 
Our approach to this problem makes use of the concept of dressed states. This formalism, originally introduced 
in [l8j] , was already employed to investigate several situations [IE, HE, HH, HE, [IE, HH Hll . It accounts for the fact 
that, for instance, a charged physical particle is always coupled to the force (gauge) field; in other words, it is always 
"dressed" by a cloud of quanta of the gauge field. In general for a system of matter particles, the idea is that the 
particles are coupled to an environment, which is usually modeled in two equivalent ways: either to represent it by 
a free field, as was done in Refs. [IE, HE], or to consider the environment as a reservoir composed of a large number 
of noninteracting harmonic oscillators (see, for instance, [28], HE, HE, HH). In both cases, exactly the same type of 
argument given above in the case of a charged particle applies, with the appropriate changes, to such systems. We 
may then speak of the "dressing" of the set of particles by the ensemble of the harmonic modes of the environment. It 
should be true in general for any system in which material particles are coupled to an environment. In atomic physics, 
the semiqualitative idea of a "dressed atom" has been largely employed in studies involving the interaction of atoms 
and electromagnetic fields (32j. In the realm of general physics, the dressing of a matter particle by an environment 
has found an application in describing the radiation damping of classical systems [33| . Our dressed states can be 
viewed as a rigorous version of these dressing procedures, in the context of the model employed here. 

We will consider our system in this paper as consisting of two atoms, each one of them interacting independently 
inside a spherical cavity with an environment provided by the harmonic modes of a field. We take it as a bipartite 
system, each subsystem consisting of one of the dressed atoms. We will consider a superposition of two kinds of states: 
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either all entities (both atoms and the field modes) are in their ground states, or just one of the atoms lies in its 
first excited state, the other one and all the field modes being in their ground states. The analysis of the (reduced) 
density matrix of the system leads to the computation of the von Neumann entropy, which measures the degree of 
entanglement of the two atoms. 

The dressing formalism for just one atom inside a cavity is briefly reviewed in Section 2 in order to establish basic 
notation and formulas for the time evolution of the states. In Section 3 the formalism is generalized for the two-atom 
system and describe the evolution of its density matrix, either in the case of a very large cavity (with infinite radius, 
that is, free space) or of a small cavity. The entanglement of the two atoms is discussed in Section 4. Finally, in 
Section 5 we present our conclusions. 



II. A SINGLE DRESSED ATOM 

Before tackling the case of two atoms, it is convenient to reproduce here the analysis of Ref. [24| for the simpler 
situation of just one atom, dressed by its interaction with the environment field. We present, in this section, a short 
review of the formalism introduced in previous works. 

We shall thus consider a atom in the harmonic approximation, coupled linearly to an environment modeled by the 
infinite set of harmonic modes of a scalar field, on the inside a spherical cavity. A nonperturbative study of the time 
evolution of such a system is implemented by means of dressed states and dressed coordinates [l8j . In particular, our 
dressed states are not the same as those currently employed in the literature, usually associated to normal coordinates. 
Our dressed states are given in terms of our dressed coordinates and allow a rigorous study of the time evolution of 
quantum systems in the context of the model employed here. The results we obtain by these means are those expected 
on physical grounds, but contain corrections with respect to the formulas obtained from perturbation theory. 

Let us start by considering an atom labeled A, having bare frequency u\, linearly coupled to a field described 
by N (— > oo) other oscillators, with frequencies u>k, k = 1, 2, . . . , N. The whole system is contained in a perfectly 
reflecting spherical cavity of radius R, the free space corresponding to the limit R — > oo. Hereafter, we shall refer to 
the harmonic oscillator as the atom, to distinguish it from the harmonic modes of the environment. Denoting by q\(t) 
(p\(t)) and qk(t) (pk(t)) the coordinates (momenta) associated with the atom and the field oscillators, respectively, 
the Hamiltonian of the system is taken as 
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where w\ is a constant and the limit N — *■ oo will be understood later on. 
principal axis by means of a point transformation, 



The Hamiltonian {T]) can be turned to 
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where //(A) = (A, {&}), k — 1, 2, . . . ,N, and r\ — 0, . . . , N , performed by an orthonormal matrix T = (t^)). The 
subscripts /i — A and /i — k refer respectively to the atom and the harmonic modes of the field and r\ refers to the 
normal modes. In terms of normal momenta and coordinates, the transformed Hamiltonian reads 
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where the Q rx 's are the normal frequencies corresponding to the collective stable oscillation modes of the coupled 
system. 

Using the coordinate transformation q^(\) — X)r A ^\\)Qr x m the equations of motion and explicitly making use of 
the normalization condition 
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with the condition 
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The right-hand side of equation ([B]) diverges in the limit N — > oo. Defining the counterterm Sto 2 — Nr] 2 , it can be 
rewritten in the form 
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Equation ([7]) has TV + 1 solutions, corresponding to the N + 1 normal collective modes. It can be shown [18| that 
if u! 2 > 5lu 2 , all possible solutions for U 2 are positive, physically meaning that the system oscillates harmonically in 
all its modes. On the other hand, when io\ < 5uj 2 , one of the solutions is negative and so no stationary configuration 
is allowed. 

Therefore, we just consider the situation in which all normal modes are harmonic, which corresponds to the first 
case above, to 2 > Suj 2 , and define the renormalized frequency 



= lim 



N- 



A - Nrii) , 

following the pioneering work of Ref. [35]. In the limit N — * oo, equation ([7]) becomes 
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We see that, in this limit, the above procedure is exactly the analogous of mass renormalization in quantum field 
theory: the addition of a counterterm —Nri\q\ (N — > oo) allows one to compensate the infinity of lo\ in such a way 
as to leave a finite, physically meaningful, renormalized frequency Q\. 
To proceed, we take the constant rj\ as 
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where Atu is the interval between two neighboring field frequencies and g is the coupling constant with dimension of 
frequency. The environment frequencies lo^ can be written in the form 
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and, so, Aw = ttc/ R. Then, using the identity 
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equation ([9]) can be written in closed form: 
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The elements of the transformation matrix, turning the atom-field system to principal axis, are obtained in terms of 
the physically meaningful quantities Cl rx and u)\ after some rather long but straightforward manipulations [18j . They 
read 
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Let us now consider the eigenstates of the system atom(A)-field, h, h, ••■}, represented by the normalized eigen- 
functions, written in terms of the normal coordinates {Q rx \i 



k xh i 2 ...(Q,t) = ll 



n. 
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where £/; s stands for the / SA -th Hermite polynomial and 
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is the normalized vacuum eigenfunction, A/a being the normalization factor. 

We introduce dressed coordinates q' x and {q^} for the dressed atom and the dressed field, respectively, defined by 
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where o> 
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{cja, o;,}. In terms of the dressed coordinates, we define for a fixed instant, t = 0, dressed states, 



\k\, «i, K2, ■ • ■) by means of the complete orthonormal set of functions [18| 
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where /z(A) labels collectively the dressed atom A and the field modes, 1, 2,3,.. ., q'^x) ~ ?A> Wi}- The ground state 
Tq in the above equation is the same as in equation p6[) . The invariance of the ground state is due to our definition of 
dressed coordinates given by equation (p~8|) . Notice that the introduction of the dressed coordinates implies, differently 
from the bare vacuum, the stability of the dressed vacuum state since, by construction, it is identical to the ground 
state of the interacting Hamiltonian ([3]). Each function ip Kxfil ...(q') describes a state in which the dressed oscillator 
q'^ is in its K^xyth excited state. 



Let us consider the particular dressed state 



^(0)) at t = 0, represented by the wave function V'oo— 1(^)0-- 



It describes the configuration in which only the /x-th dressed oscillator is in the first excited level, all other being in 



their ground states. As shown in Ref. [lj|, the time evolution of the state 

r5 t(A) (t))=E/^(*)| r i (A) (°))» 

v 

where //(A), f(A) = A, {z}, with {i} referring to the field modes, and 
Moreover, it can be shown that, for all ^, 

Ei/^(t)i 2 = i, 



is given by 



(20) 
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which allows to interpret the coefficients f^{t) as probability amplitudes; for example, f\\(t) is the probability 
amplitude that, if the dressed atom is in the first excited state at t = 0, it remains excited at time t, while fxi(t) 
represents the probability amplitude that the i-th dressed harmonic mode of the field be at the first excited level. 



III. TIME EVOLUTION OF A DRESSED TWO-ATOM STATE 



We now consider a bipartite system composed of two subsystems, A and 6; the subsystems consist respectively of 
dressed atoms A and B, in the sense defined in the preceding section, the whole system being contained in a perfectly 
reflecting sphere of radius R. Let us consider the eigenstates of the subsystems A and B with X — A, B labeling the 
quantities referring to the subsystems. 
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We consider the Hilbert space spanned by the dressed Fock-like states, 



,}j = \n A ,k 1 ,k 2 , n B ,qi,q 2 , ...)= \T$ A)ku k 2 ,..) 



(23) 



in which the dressed atom A is at the tia excited level and the atom B is at the n B excited level; the (doubled) 
dressed modes of the field are at the ki, k%, ■ ■ ., qi, <fe, • • ■ excited levels. Using this definition, let us consider at time 
t = 0, a family of entangled states of the bipartite system given by 
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where < £ < 1. In equation ([24l 



r i(A)o(B)oo---( )) and 
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the dressed atom A (B) is at the first level, the dressed atom B (A) and all the field modes being in the ground state. 
They are 
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The density matrix at t = is 



e (o) = |¥>(*| 

= £ |U, 0, 0, ■ • ■ ; B , 0, 0, • ■ •} (1 A , 0, 0, • • • ; Q B , 0, 0, • • -| 

+ (1 - |0 A , 0, 0, • • ■ ; 1 B , 0, 0, ■ • ■} {0 A , 0, 0, • • • ; 1 B , 0, 0, • ■ -| 
+ I U, 0, 0, • • ■ ; B , 0, 0, • • ■) (0 A , 0, 0, • • ■ ; 1 B , 0, 0, ■ 

+ ^£(1-0^ |0a, 0, 0, • • ■ ; 1b, 0, 0, • ■ •) (U, 0, 0, • • • ; B , 0, 0, ■ • 
At time t, the state of the system is described by the density matrix 

g(t) = er lHt \V) (*|e 
where -ff is the Hamiltonian of the whole system, such that 
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and Ha and H B are the Hamiltonian H\ of equations (fTJ) or ([3]) . We then obtain 
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where the states |r^ 00 ...), |r^ 00 ...) are stationary and the states |r^ 00 . ..(£)), |r^ ...(i)) evolve according to equation 
(®. 

In order to investigate how the superposed states evolve in time, we shall consider the reduced density matrix 
obtained by tracing over all the degrees of freedom associated with the field. The computation is analogous to the 
one presented in Ref. [24| . After taking the trace, the density matrix has the indices referring to the 2-atom states. 
Explicitly, we have 
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In the above expression we have typically 

(n A , fci, fc 2 , . . . |r^ 00 ...) = 4!o4 2 o • ■ 
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so that the sums in the elements of the reduced density matrix are of one of the types below: 
E (nA,k l7 k 2 , . . . |r^ 00 ...) (ro 00 ...|m^,fci,fc 2 , • • •) = 
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Collecting all these expressions and their analogues for the atom B into the elements of the reduced density matrix 
we finally obtain 
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That is, the nonvanishing elements are given by 



PoTilW = (W)IW*)I 



(37) 



PifoK*) = V^T)e^fAA(t)r BB (t), 

where equation (122p was used. We check immediately that the trace of this reduced density matrix is one, 

P&Z+pZll+PllZV+Pullit) = !• (38) 

This property ensures that p represents physical states of the system. Also, we see that Tr[p 2 ] ^ 1 and therefore, the 
superposed states are not pure. The degree of impurity of a quantum state can be quantified by the departure from 
the idempotency property. In the present case: 



£>(t,0 = 1-Tr[p 2 ] 



2 (e \fAA(t)\ 2 + (1 - \fBB(t)\ 2 ) ~ 2 (e \fAA(tf + (1 - \fBB{t)\ 



(39) 



In the remainder of this section we consider the two atoms as identical and, accordingly, we adopt the subscript 
for both of them. A = A = B = 0; we also define 



9A = 9b = g; rjA = riB = r l> ua = wb = w', /aa(*) = /bbW = /oo(*)- 



(40) 



In this case, the matrix elements in equations (|37[) simplify and, from equation (|39p . we see that the degree of impurity 
becomes independent of the superposition parameter £: 



D(t,0=2\f 00 (t)\ 2 (l-\f o(t)\ 2 ). 



(41) 



In order to pursue the study of the time evolution of the superposition of the two-atom states, we have to determine 
the behavior of /oo (0- We shall analyze it in the situations of a very large cavity (free space) and of a small one. 



A. The limit of an arbitrarily large cavity 



We start from the matrix element t 7 ^ in equation (| 14[) and consider an arbitrarily large radius R for the cavity. 
The two atoms behave independently from each other, so let us focus on just one of them, either the atom A or the 
atom B. Remembering that i] = y/ Age/ R, we have 



Urn tg= lim _^M^B2=. (42) 
R ^°° R ^°° ^(ft 2 - Co 2 ) 2 + 4g 2 n 2 

In this limit, Alu = nc/R —> cLu = dfl and the sum in the definition of /oo(i), equation (|21| . becomes an integral, so 
that 

/oo(*) = ^/ dSl 6 (43) 

n Jo (n 2 - to 2 ) + Ag 2 n 2 
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We then proceed as in [24|]. We define a parameter k = y/u; 2 — g 2 and consider whether k 2 > or k 2 < , for 
which n 2 ^> and k 2 <0 correspond respectively to weak (j < (ii) and strong [g u>a) coupling of the atoms with 
the environmment. For definiteness we consider in the following the weak-coupling regime. We get in this case [HI 



-gt 



/oo(*) = e 

where the function G(t;uj,g) is given by 

G(t;Q,g) 
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For large times, the quantity \foo(t)\ 2 is given by 2 
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As t — » oo, we see that the expression for |/oo(^)| 2 §° to zero. 



B. Small cavity 



For a finite (small) cavity, the spectrum of eigenfrequencies is discrete, Alu is large, and so the approximation 
made in the case of large cavity does not apply; no analytical result can be obtained for /oo(i) in this case. For a 
sufficiently small cavity, the frequencies Q r can be determined by following the steps described in [24j. Let us label 
the eigenfrequencies as flo, {fife}, k — 1,2,... Then, defining the dimensionless parameter 



.9 = gR 

AlU 7TC ' 
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we rewrite equation (|13p in the form 
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Taking 6 <C 1, which corresponds to R <C 7rc/(7 (a small cavity), it is shown in [24[ that, for k = 1, 2, . . ., the solutions 
are 
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If we further impose that 5 < 2g 2 /ttlu 2 , a condition compatible with (5 <C 1, then fig is found to be very close to lu, 
that is, 



fio ~ w 1 



7T<5 



(50) 



To determine foo(t), we have to calculate the square of the matrix elements and (tfy . They are given, to 

first order in S, by 
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We thus obtain, for sufficiently small cavities (S <C 1), 
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To order 6 2 , a lower bound for |/oo(^)| 2 is obtained by taking the value 
using the tabulated value of the Riemann zeta function £(2) = 7r 2 /6, 



T for both cosines in the above formula, 



|/oo(*)| 2 > [l + ^S 



4tt5 
~3~ 



4tt 2 5 2 
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We see that the quantity |/oo(i)| 2 , which dictates the behavior of the density matrix elements and of the measure 
of purity in equation (|41[) . has very different behaviors for free space or for a small cavity. This implies that in the 
situation of a small cavity, in contrast to the free space case, all matrix elements in equations (|3T[) are different from 
zero for all times. 

In Figure (|TJ) the degree of impurity from equation (|4"Tj) is plotted as a function of time in the cases of an arbitrarily 
large cavity (R — > oo) and of a small cavity. We take S — 0.1, with Q — 1.0 and g = 0.5 fixed (in arbitrary units). 




FIG. 1: Behavior of the degree of impurity D as function of time, equation (|4ip . for a small cavity (dashed line) and a very 
large cavity (solid line); we take the parameters g = 0.5, 8 = 0.1 and ui = 1.0 (in arbitrary units). 

We see from the figure that for a very large cavity (free space) the two-atom system dissipates; with the passing of 
time, both atoms go to their ground states. For a small cavity the system never completely decay. 



IV. TIME EVOLUTION OF THE ENTANGLEMENT 



In order to study how the entanglement of the two-atom states evolves in time, we shall, as before, consider the 
system as a bipartite system, in which each atom carries its own dressing field. In this way, we shall determine the 
time behavior of the von Neumann entropy associated with the reduced density matrix with respect to one of the 
subsystems, which is given by taking the trace over the states of the complementary subsystem. 

Let us initially treat the system at t = 0. Then we have, for the subsystem A, the reduced density matrix 

Pa(0) = TMI*>(*I) 

OO 

= C|1a,0,0,---)<1 A) 0,0,---| + (1-O|0a,0,0,---)<0 a ,0,0,---| (54) 
and, similarly, for the subsystem B, 

pb(o) = Tr A mm 

= C |0 B ,0,0, •• •}(0 B ,0,0,---| + (1 -Oils, 0,0, ••■}(1b,0,0,---|. (55) 

The degree of entanglement of the two-atom system is measured by the von Neumann entropy of any of the reduced 
density matrices; for instance, 

E(0 = -Tr [p A In p A ] = - £ a In a, (56) 

a 

where the sum is taken over the eigenvalues a oi p A . Since p A is diagonal in the Fock basis of the dressed states of 
the atom A, its eigenvalues can be read directly from ([54]) : 



«i = 1 — £, a 2 — £, tti = Q!3 = ■ • • = 0. 



(57) 
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Therefore, 



£(0 = -[(i -0 Mi -0 + 



(58) 



The time evolution of the states 



MB) 
1(A)0(-B)00- 



and 



(AB) 

0(A)1(B)00- 



are governed by the time evolution of the states 



|r^ 00 ...) and |rf 00 ...), respectively, given by equation (f2"0"|) . 

|r?oo--(*)) = S)/Av(*) 



(59) 



where, in accord with the notation of the preceding section, the label A now refers to each one of the dressed atoms 
A and B and 



-itt.t 



(60) 



In equation (|59[) . -Tioo— ^ i s the state in which the dressed mode v{\) of the atom A is at the first level and all the 

other dressed modes are in the ground state. 

The reduced density matrix corresponding to the subsystem A at time t is 

PA (t) = Tir B Q(t)=Tr B [\*(t))(*(t)\}. 

Using equation (|5"9"|) , one writes p^(t) in terms of the quantities f\v{t) from equation (f6"0|) : 



PA(t) = E ( n B;1i><l2, ■■■ |*(*)> (*(*)! n B ,q 1 ,q 2 , .. .) 



E E [Ve/A^w 



100- 



{n B ,qi,q2, ■ ■ ■ |rf 00 ...) + \j\ - ie^fB^t) |r^ 00 ...) (n B ,qi,q 2 , 



MB) 
100- 



X 



[VftXX*) (r$£| <r s 00 ...| n B , gi ,g 2 , . . .) + A /We- i ^ I/ (t) (i 

££M*)/JU*) | r io^-) ( r ioo. } .| + (i - |r&o.„> <r^ 



"6,91,92, ■ 



000--- ' 



where we have used 



r i0^ } - ) = fBB{t)S nB i Y\ S i0 + 5 nB Q ^ /si(*)<**l H <*j'0 



(61) 



(62) 



and equation 

The time-dependent von Neumann entropy is now given by 



E(t,0 = -Tr \p A (t) \n PA (t)\ = -E aln «> 



(63) 



where here a are the time-dependent eigenvalues of the reduced density matrix. These should be solutions of the 
so-called characteristic equation, which in the case of l|6ip. reads 



det 



£|/aa| ~Oi UA2f*AA 
ifAAf'* Al t\fAl\ 2 -a if A 2f* A1 

£/aa/A 2 £/ai/I 2 Z\fA2\ 2 -a 



V 



We thus find that the nonzero eigenvalues of p A are 



(64) 



(65) 
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This then implies that the von Neumann entropy takes the expression 

E(t, = - [(1 - HI - + € In (0] , (66) 

that is, all the time dependence of the von Neumann entropy for this two-atom system, coming from the f\ v (t), 
is completely cancelled in the computation of the entropy, in all situations, thereby reproducing exactly the same 
expression as in the t = case, with the maximum entanglement occuring at £ = 1/2 (see Figure [2]). In other words, 
although the superposition of states evolves in time, in different ways in the limits of a very large cavity and of a 
small one, the entangled nature of these two-atom states remains unchanged for all times, independently of the size 
of the cavity. 

E 

0.7 
0.6 
0.5 
0.4 
0.5 
0.2 
1 

0.2 0.4 0.6 0.8 1.0 3 

FIG. 2: Behaviour of the von Neumann entropy at all times, equation (|66[) . as a function of the parameter £. 



V. CONCLUDING REMARKS 

In this paper we have considered a system composed of two atoms in a spherical cavity, each of them in independent 
interaction with an environment field. The model employed is of a bipartite system, in which each subsystem consists 
of one of the atoms dressed by the field. We make the assumption that initially we have a state in which one of 
the dressed atoms is in its first excited level and the other atom and the field modes are all in the ground state, is 
superposed with a state in which the atoms have their roles reversed. 

The time evolution of the superposed states leads to a time-dependent (reduced) density matrix. Expressions for its 
elements are provided in both the cases of an infinitely large cavity (that is, free space) and of a small one, when the 
two atoms are considered as identical. Very different behaviors are obtained for this time evolution. In the large-cavity 
case, the system shows dissipation, and, with the passing of time, both atoms go to their ground states. For a small 
cavity, an oscillating behavior is present, so that the atoms never fully decay. 

Nevertheless, in spite of these rather contrasting behaviors and of the nontrivial time dependence of the density 
matrix, we obtain a time-independent von Neumann entropy, which means that the initial entanglement of the two 
atoms remains unchanged as the system evolves. 
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